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We study dynamics of a phase boundary in a one-dimensional lattice gas, which is initially put 
into a non-equilibrium configuration and then evolves in time by particles performing nearest- 
neighbor random walks constrained by hard-core interactions. Initial non-equilibrium configu- 
ration is characterized by an 5"-shape density profile, such that particles density from one side 
of the origin (sites X < 0) is larger (high density phase, HDP) than that from the other side 
(low-density phase, LDP). We suppose that all lattice gas particles, except for the rightmost 

^ particle of the HDP, have symmetric hopping probabilities. The rightmost particle of the HDP, 

which determines the position of the phase separating boundary, is subject to a constant force 
F, oriented towards the HDP; in our model this force mimics an effective tension of the phase 

■i ■ separating boundary. We find that, in the general case, the mean displacement X(t) of the phase 

boundary grows with time as X(t) = a(F)t 1 ' 2 , where the prefactor a(F) depends on F and on 
the initial densities in the HDP and LDP. We show that a(F) can be positive or negative, which 
means that depending on the physical conditions the HDP may expand or get compressed. In 
the particular case when a(F) = 0, i.e. when the HDP and LDP coexist with each other, the 
second moment of the phase boundary displacement is shown to grow with time sublinearly, 



C3 



X 2 (t) = jt 1 ' 2 , where the prefactor 7 is also calculated explicitly. Our analytical predictions are 



shown to be in a very good agreement with the results of Monte Carlo simulations. 

o 
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I. INTRODUCTION. 

o 

A fundamental question concerning the behavior of systems out of equilibrium can be 
formulated as follows. Suppose that two different phases, composed of the same or of 
two different substances, are initially prepared in different regions of space and have a 
common interface. What is the future evolution for the system and for the phase-separating 
interface? This problem appears in the analysis of such diverse phenomena as expansion of 
the poisoned state in catalytic reactions, or, more generally, propagation of chemical fronts, 
dielectric breakdown, growth of dendrites and clusters in Ising magnets, spatial intermittency 
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in hydrodynamics, wetting, rise of liquids in capillaries and many others (see [p]-Tl^| and 
references therein). 

Theoretical analysis of the problem follows basically two distinct avenues. One type of 
approach is to describe the system evolution in terms of some appropriate set of starting 
equations, a standard list of which includes such non-linear differential field equations of 
different complexity as, e.g. Newell- Whitehead, viscous Burgers, Swift-Hohenberg, Cahn- 
Hilliard equations and etc |PHT^|. These equations are usually referred to as "microscopic", 



with the understanding, however, that they don't involve atomic degrees of freedom but 
rather serve as elementary building blocks from which the analysis starts. Another ap- 
proach consists in the direct study of models involving particles with microscopically defined 
dynamics |T3|-[16|. The dynamic rules can be, for instance, chosen to construct a suitable cel- 



lular automaton, which converges asymptotically to the field equation in question, e.g. the 



Navier-Stokes equation fl5|-|18f, allowing then for much more efficient numerical analysis than 
simulations of the continuous-space counterpart. Alternatively, they can be deduced from 
realistic microscopic interactions with the intention to derive equations describing the time 



evolution of some macroscopic properties, such as, e.g. local particle densities [0,|T3-pT 



Considerable progress has been made recently in this direction [jT5| , pT5| ^pT[] , which revealed, 
however, the fact that macroscopic equations derived on the basis of realistic microscopic in- 
teractions may have a different structure compared to the generally accepted field equations 
and can be reduced to them only under certain assumptions. 

In the present paperf] we study dynamics of the phase boundary in a two-phase micro- 
scopic model system consisting of identical hard-core particles which are placed initially 
on a one-dimensional, infinite in both directions lattice in a non-equilibrium, "shock" -like 



*This paper is based partly on the talk given at the conference on Inhomogeneous Random Sys- 
tems, Palaiseau, France, January 1997, and at the workshop on Instabilities and Non-Equilibrium 
Structures, Santiago, Chili, December 1997 



configuration. That is, particles mean densities from the left and from the right of the 
origin of the lattice (Fig.l), which we denote as p_ and p + respectively, are generally not 
equal to each other. We suppose, without lack of generality that p_ > p+ > 0, and will 
call in what follows the phase which initially occupied the left half-line as the high-density 
phase (HDP), while the phase initially occupying the right half-line will be referred to as the 
low-density phase or the LDP. Particles are then allowed to perform symmetric, (i.e. with 
equal probabilities for going to the left or to the right), hopping motion between the nearest 
lattice sites under the constraint that neither two particles can simultaneously occupy the 
same lattice site and can not pass through each other. Further on, we single out the right- 
most particle of the HDP, which determines position of the phase boundary and which we 
will call in the following as the PBP - the "phase boundary" particle (Fig.l). We suppose 
that this only particle is subject to a constant force F which favors the PBP to jump in a 
preferential direction. Thus for the PBP the probabilities of going to the right (p) and to 
the left (q) will be different from each other. In most situations we will suppose that F is 
directed towards the HDP; we adopt the convention that in this case F is positive definite, 
F > 0; the PBP hopping probabilities p and q are related to the force and the reciprocal 
temperature j3 through p/q = exp(—f3F) and p + q = 1. From the physical point of view, 
such a constant force can be understood as an effective boundary tension derived from the 
solid-on-solid-model Hamiltonian of the phase-separating boundary |25-27| and mimic, in a 



mean-field fashion, the presence of attractive interactions between the lattice-gas particles 
which are not explicitly included into the model. 

We hasten to remark that the system evolution in the case when long-range attractive 
interactions between the gas particles are present can be fairly more complex. First, in this 
case the hopping probabilities of any given particle are coupled to the instantaneous positions 
of all other gas particles and, consequently, evolution of the local densities is described by 
non- linear and non-local integro-differential equations (see ||T9| , |20|| and references therein). 



Furthermore, the "boundary-tension" force F is generally time-dependent and reaches a 
constant value only at sufficiently large times. Moreover, this value will be explicitly depen- 



dent on the density distribution around the PBP resulting in a non-linear coupling between 
the dynamics of the PBP and of the lattice-gas particles. A mean-field-type analysis of 
this situation, which has interesting applications within the context of spreading of liquid 
monolayers on solid supports [p2] , P3|j , has been presented recently in P5| ; here, we assume 
that F is a fixed given parameter, which is independent of the particle density. We note 
that such an assumption is justified when particle-particle interactions are sufficiently weak 
and time t is sufficiently large [^ . 

We note now that several particular cases of the general model under study have been 
already discussed in the literature. Consequently, considering different limiting with respect 
to p + , p_ and F situations we will be able to check our predictions against already known 
results. We mention some of these results: 

(a) Dynamics of the PBP in the symmetric case p_ = p + = p and in the absence of an 
external force F has been studied as early as 1965 by Harris |24j| , who has shown rigorously 



that the mean-square displacement of the PBP grows sublinearly with time, 

i-p \n 



Here and henceforth the overline denotes the averaging with respect to different realizations 
of the PBP trajectories X r {t). Moreover, it was shown in |29|] that t^ l ^X r {t) converges in 



distribution to a Gaussian variable with variance (1 — p)J2/n/p. 

(b) Further on, a rigorous probabilistic description of the situation with p_ = 1, p + = 



and zero "boundary tension" force has been developed in ||29|| . It was proven that in this 

case the mean displacement of the PBP grows in time in proportion to Jt log{t) as t — > oo. 

(c) An opposite case when F = oo, (such that the PBP performs totally directed random 

walk), while p + = and p_ = p has been considered in |J{J. It was demonstrated that the 



mean displacement of the PBP obeys X r {t) = —au m \/t, in which law the prefactor a>u m is 
defined implicitly by 

^1 a Um exp(«L/2) [1 - $(a Um /y/2)] = 1 - p, (2) 



where <&(x) denotes the error function. 



(d) A more general situation has been considered in [[!![] and subsequently, in ||32| , which 
works deal with the behavior of the driven PBP in the symmetric case p_ = p + = p. Both 
works have shown that at arbitrary negative values of the boundary tension force, F < 0, 
the mean-square displacement of the SFB follows 

X^T) = a(F) Vi (3) 

In [[J]]] it was found analytically, in terms of a mean-field-type approach, and also confirmed 
by numerical Monte Carlo simulations that for arbitrary p and p > q the parameter a(F) is 
determined by the following transcendental equation: 

(yf a(F) exp(a 2 (F)/2) [1 + *(a(F)/>/2)] + 



+ y - HK HJ )Urk «(*") ex P (« 2 (F)/2) x 



p -qiX-p) in 
p-q K V2 



:[1 - *WF)/V2)] + iZlfLlA) = ^L ( 4 ) 

p — q [P — Q) 



In |j32|1 ) which has also established an interesting relation between the time evolution of a 
symmetric lattice gas with a single driven tracer and the evolution of the interfaces in a two- 
dimensional Potts model with Glauber dynamics, the PBP dynamics was analysed in terms 
of a rigorous probabilistic approach and the result in Eq.(4) has been rigorously proven. 
An interesting observation made in |5T| and subsequently, in [52|, concerned the validity 



of the Einstein relation for the tracer diffusion in a one- dimensional hard-core lattice gas. An 
illuminating discussion of this issue and a considerable amount of new results establishing 
the Einstein relation for different interacting particle systems was presented recently in J53 . 



Now, Eq.(4) shows that in the case of a weak asymmetry, i.e. when p — q — > 0, the parameter 
a(F) is given exactly by 

a(F) = (p - q) l -^-P- J| (5) 
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If one defines then the time-dependent mobility //(£) of the driven PBP as //(£) = 
lirriF-+oX(t) / Ft, it would yield, by virtue of Eq.(3), that /j,(t) = t^ 1 ^ 2 limF_^ a(F) / F , where 
ct(F) is given by Eq.(5). On the other hand, the time-dependent diffusivity D(t) of the PBP 
in absence of external force obeys D(t) = (1 — p)/p(27it) 1 ^ 2 ||24| . As it was shown in f31 



and [P^ ], the asymptotic form in Eq.(5) implies that /j,(t) = (3D(t), i.e. the Einstein relation 
holds exactly in the non-stationary regime^]. 

Here we develop an analytical dynamical description of the lattice gas model with initial 
"shock"-like configuration of particles, which is based on the mean-field-type assumption 
that the average of the product of local realization-dependent variables, describing occupa- 
tion of lattice sites, factorizes into the product of their average values, i.e. the local densities. 
Such an assumption permits us to derive the closed-form system of equations for the time 
evolution of the PBP mean displacement and of the density profiles around it. These equa- 
tions allow for the analytical solution and explicit computation of the time evolution of the 
PBP mean displacement. Our main results are the following: 

We find that at sufficiently large times the mean displacement of the PBP obeys Eq.(3), 
in which the prefactor a(F) is determined implicitly as the solution of the transcedental 
equation 



qp-{l + Jn/2a(F) exp(a 2 (F)/2) [1 + sgn(a(F)) ^(^(F)!/^)]}- 1 



-pp + {1 - s/n/2a(F) exp(a 2 (F)/2) [1 - sgn(a(F)) ^HF)!/^)]}" 1 = q - p, (6) 

where sgn(a) = 1 for a > and sgn(a) = —1 for a < 0. Eq.(6) holds for any rela- 
tion between p, q and p± and reduces to the previously obtained Eq.(2) and Eq.(4) in the 
appropriate limits. This result has been found subsequently in |32|. 



^Two of us, G.O. and S.F.B., wish to thank Professor J.L.Lebowitz who has suggested us to 
examine the question of validity of the Einstein relation in the non-stationary regime. 



Equation (6) predicts that three different regimes can take place depending on the rela- 
tion between p/q and p±: 

(1) When p(l — p + ) > q(l — p_) the parameter a(F) is finite and positive definite, which 
means that the HDP expands compressing the LDP. In the particular case p/q = 1 and 
p + = the parameter oe(F) appears to be a positive, logarithmically growing with time 
function, which behavior agrees with the results of |29| . 



(2) When p(l — p + ) = q(l — p_), the parameter a(F) = 0. This relation between the 
system parameters when the HDP and the LDP are in equilibrium with each other and the 
PBP mean displacement is zero, was found also in J33| from the analysis of the stationary 
behavior in a finite one- dimensional lattice gas (see also Section VI. A). Despite the fact that 
the PBP mean displacement is zero, the fluctuations in the PBP position grow with time. 
We show here that in this case the mean-square displacement of the PBP obeys 

(l-p_)(l-p + ) [to 



(P- + P+ - 2 p_ p+) V vr 
which reduces to the classical result in Eq.(l) in the limit p = q and p_ = p + . Eq.(7) is 
derived here using heuristic arguments based on the Einstein relation between the mobility 
and diffusivity of a test particle and is confirmed by Monte Carlo simulations. 

(3) When p(l — p + ) < q(l — p_) the parameter a(F) is less than zero - the expanding 
LDP and the applied force effectively compress the HDP. 

Further on, we show that the particles density profile as seen from the moving PBP, 
stabilizes around its position, approaching a constant value. We find that in the regime 
when the HDP expands (i.e. a(F) > 0) the density profile p(X; t) around the PBP is given 
by 

p(X < X(t); t) « 1 + /+ P (KF)|) [1 + a(F) r x ' 2 (X{t) - X) + ... ], (8.o) 

p{X > X(t);t) « — ^j-— -- [1 + a(F) r 1 / 2 (X (t) - X) + ... ], (8.6) 

for X(t)—X <^ y/i/A. Within the opposite limit {cx(F) < 0) when the HDP gets compressed, 
p(X; t) follows 



p(X < X(t); t) » - /" ,„ [1 + a(F) t~^(X(t) - X) + ... ], (9.o) 
1- J_(|a(F)|) 



p(X > X(t);t) » 1 + /+ P ( | ft(F)|) [1 + o(F) t-V 2 (X(t)-X) + ... ], (9.6) 



where 



a 2 (F) 



/±(|a(F)|) = a\F) rdzexp(-^-^(z 2 T 2z) = 

Jo 2 

= y||a(F)| exp(^) [1 ± ma(F)\/V2)] (10) 

This paper is outlined as follows: in Section 2 we formulate our model and discuss the 
approximations involved. In Section 3 we write down basic equations, describing the time 
evolution of the PBP and of the lattice-gas particles in the particular case when the LDP 
is absent. Section 4 presents solutions of the dynamic equations in the case p + = and 
discussion of the PBP dynamics at different values of the initial mean density p_ and of the 
"boundary tension" force F. Further on, in Section 5 we consider the general case when the 
LDP is present and evaluate the dynamical equations describing the time behavior of the 
system under study. Section 6 is devoted to the analysis of these equations and evaluation 
of their solutions. Next, Section 7 contains the description of the Monte Carlo simulations 
algorithm. Finally, in Section 8 we conclude with a brief summary of results and discussion. 

II. THE MODEL. 

Consider a one-dimensional, infinite in both directions regular lattice of unit spacing, the 
sites {X} of which are partly occupied by identical particles. Suppose next that the initial 
configuration of particles is as depicted in Fig.l, i.e. all particles are placed on the lattice 
with the single occupancy condition, at random positions and in such a way that the mean 
particle density (p) at sites X < is different from the mean particle density (p + ) at the 
sites with X > 0. As we have already mentioned, we supposed that p_ > p + > and call 



the gas phase from the left of the PBP as the HDP and the gas phase from the right of the 
PBP as the LDP. 

After deposition onto the lattice, the particles are allowed to move by attempting jumps 
to neighboring sites. The motion is constrained by the hard-core exclusion between the 
particles; that is, neither two particles can simultaneously occupy the same lattice site nor 
can pass through each other. For a given realization of the process the instantaneous particle 
configuration is described by an infinite set of time-dependent occupation variables {rx(t)}, 
where each Tx(t) assumes two possible values; namely, Tx(t) = 1 if the site X is occupied 
at time t and Tx(t) = if this site is vacant. Position of the PBP at time t is denoted as 
X r (t), which is also a random, realization-dependent function. 

More specifically, we define particles dynamics as follows: each particle waits a random, 
exponentially distributed time with mean 1 and then selects, with given probabilities, the 
direction of jump - to the right or to the left. When the direction of the jump is chosen, 
the particle attempts to jump onto the nearest site. If the target site is unoccupied by any 
other particle at this moment of time, the jump is instantaneously fulfilled. If the site is 
occupied, the particle remains at its position and waits till the next attempt. The process 
is memory-less and the choice of jump directions for different attempts is uncorrelated. 

We will distinguish between the jump probabilities of the PBP and the jump probabilities 
of all other particles of the gas. The jump probabilities of the gas particles are symmetric, 
i.e. for them an attempt to jump to the right and an attempt to jump to the left occur with 
probability 1/2, while the jump probabilities of the PBP are asymmetric; it attempts to 
jump to the right with probability p and to the left with probability q, p + q — 1. These 
probabilities are related to the "boundary tension" force and the temperature T = 1/(3 
through the relation p/q = exp(—[3F). 

Here we will develop a mean-field-type description of the time evolution of the system 
under study using an approximate picture, based on two assumptions: 

We assume first that the average of the product of the occupation variables Tx(t) of 
different sites factorizes into the product of their average values, which corresponds to the 
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local equilibrium assumption^. Under such an assumption we can describe the system evo- 
lution directly in terms of the realization-averaged values p(X; t) =< Tx(t) >, which define 
the local density of the site X at time t or, in other words, the probability that the site X 
is occupied by a lattice gas particle at time t. The resulting equations will then be closed 
with respect to p(X;t), i.e. will not include higher-order correlation function. A non-trivial 
aspect of these equations, which is common for diverse front propagation problems @, is 
that one of the boundary conditions is imposed in the moving frame. 

Secondly, the evolution of the spatial position of the phase boundary will be described 
in terms of P(X;t), which defines the probability of having the PBP at site X at time t. 
Anticipating that in the limit t — > oo the ratio 



XJt) , s 

r H= -> oo, 11 



'xm 

i.e. that fluctuations in the PBP trajectory grow at essentially slower rate than its mean dis- 
placement, we will neglect fluctuations in the PBP trajectories, supposing that the position 
of the PBP at time t is a well-defined function of time, which is the same for all realizations 
of the process, X r (t) = X(t). We note that such an assumption is quite consistent with 
rigorous results presented in |32|]. 

These two simplifying assumptions will allow us to determine explicitly the dynamics 
of the PBP and to calculate the density profiles as seen from the PBP. Our analytical 
predictions will be checked against the results of Monte Carlo simulations of the stochastic 
process, described in the beginning of this Section, which allows a direct control of the 
validity of the local equilibrium assumption. 



^Note, that the density profiles around the PBP, as shown by Eqs.(8) and (9), tend to constant 
values on progressively larger and larger scales as the PBP advances; that is, there is no structure 
in the density profiles and they are merely the product measures. This "propagation of local 



equilibrium" [34] insures that the decoupling procedure involved is correct in a large time scale. 
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III. DYNAMICAL EQUATIONS IN THE ABSENCE OF THE LOW-DENSITY 

PHASE. 

We consider first dynamics of the PBP in the particular case when the LDP is absent 
(Fig. 2), i.e. p + = 0, and thus the jumps of the PBP away from the HDP are unconstrained. 

We start with the derivation of the dynamical equations, which govern the time evolution 
of p(X;t). In the continuous-time limit and under the assumption of the factorization of 
the occupation variables at different sites, dynamics of p(X; t) is guided by the following 
balance equation 

p(X;t) = -l p (X;t){l - p(X + l;t) + 1 - p(X-l;t)} + 

+ \ (1 " P(X; t)) {p(X + 1; t) + p(X - 1; £)}, (12) 

where the terms in the first two lines describe the contribution due to the jumps from the 
occupied site X onto the neighboring unoccupied sites, while the terms in the third line 
account for possible arrivals of particles to unoccupied site X from the occupied adjacent 
sites. One may readily notice that in Eq.(12) non-linear terms cancel each other and it 
reduces to the discrete-space diffusion equation 

p(X;t) = I{ p (X + l;t) + p(X-l;t)- 2 p(X;t)} (13) 

We note now that in deriving Eq.(13) we have implicitly supposed that the jump prob- 
abilities of particles arriving from the site X + 1 are symmetric. This means that Eq.(13) 
holds only for the sites X, which are inaccessible for the PBP, whose jumping probabilities 
are asymmetric by definition. Consequently, Eq.(13) is valid only for the sites X such that 
X < X(t) - 1. For p(X(t) - 1; t) we have instead of Eq.(12) 

p(X(t)-l;t) = l -{p{X{t)-2-t){l - p(X(t)-l;t))- 

-p(X(t)-l;t)(l - p(X(t)-2;t))}+ 
11 



+ q (1 - p(X(t) - 1; t)) p{X{t) - 2; t) - p p(X(t) - 1; t) (14) 

The terms in the first two lines of Eq.(14) describe exchanges of particles between the sites 
X(t) — 1 and X(t) — 2. The particles which may be involved in these exchanges are the 
lattice gas particles which have symmetric jumping probabilities and here, again, the non- 
linear terms cancel each other. Two terms in the third line of Eq.(14) account for the 
effective change in the occupation of the (X(t) — l)-site due to jumps of the PBP and are 
defined in the frame of reference moving with the phase boundary. The first term describes 
creation of a vacancy at a previously occupied site (X(t) — 1) due to the unconstrained jump 
of the PBP away of the gas phase. The second one accounts for the effective creation of a 
particle at (X(t) — 1) in the event when the PBP jumps onto unoccupied site (X(t) — 1) 
and the site (X(t) — 2) is occupied prior to the jump. 

Similar reasonings yield the following equation for the time evolution of the probability 
distribution P(X;t), which defines the PBP dynamics, 

P(X; t) = - P(X; t) { p + q (1 - p(X(t) - 1; t))} + 

+ p P(X - 1; t) + q (1 - p(X(t); £)) P(X + 1; t) (15) 

Multiplying both sides of Eq.(15) by X and summing over all lattice sites we find that the 
displacement of the PBP obeys 

X(t) = p - q + g/(l;*), (16) 

where we took into account the normalization condition J2x P(X', t) — 1 and denoted as 
f(X;t) the pair-wise correlation function 

f(X;t) = J2P(X;t)p(X-X;t) (17) 
x 

Equation (17) defines the probability of having at time moment t a particle at distance A 

from the PBP, or, in other words, can be interpreted as the density profile as seen from the 

moving PBP. 
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We now turn to the time evolution of f(X;t). Differentiating the pair-wise correlation 
function in Eq.(17) with respect to time, we have 

/(A; t) = £{p(X - A; t) P(X; t) + P(X; t) p(X - A; t)} (f 8) 

x 

We notice now that again the behavior for A = 1 and A > 1 has to be considered separately. 
In the domain A > 1 we find, taking advantage of Eqs.(13) and (15), that f(X;t) obeys 

/(A;t) = I{/(A-l;t) + /(A + l;i) - 2/(A;t)}- 
-( p - q )f(X;t) + pf(X-l;t) + qf(\ + l;t)- 

-q J2 P (X;t)p(X-l;t)p(X-X-l;t) + 

x 

+ q^2P(X;t)p(X-l;t)p(X-X;t) (19) 
x 

We proceed further on making the same simplifying assumption, which underlies the deriva- 
tion of Eqs.(12) to (15), i.e. assuming that the average of the product of the occupation 
variables decouples into the product of the average values. This means that the two last 
terms in Eq.(19) can be rewritten as 

J2P(X;t)p(X-l;t)p(X-X-l;t) = 

x 

= {]T P(X; t) p(X - 1; *)} {£ P(X'; t) p(X' - X - 1; *)}, (20.a) 

X X' 

and 

Y,P(X;t)p(X-l;t)p(X-X;t) = 

x 

= {]T P(X; t) P (X - 1; t)} (E P(X>; t) p(X> - A; t)} (20.6) 

X X' 

Decoupling of the third-order correlation functions as in Eqs.(20) permits us to cast Eq.(19) 
into the following form 
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/(A;0 = I{/(A-l;t) + /(A + l;i) - 2/(A;t)}- 
-/(A;t) + p/(A-l;t) + g /(A + l;f)- 

-qf(l;t)(f(X + l;t) - /(A;0) (21) 

Equation (21) does not include now the third-order correlation functions and thus is closed 
with respect to f(X;t). 

Next, using Eqs.(18), (15) and (14) and decomposing the third-order correlation functions 
into the product of pair- wise correlations, we obtain for the time evolution of /(l; t): 

/(i;*) = \if{W) - /(i;*)}- 

-f(l;t) + 2qf(2;t)(l - f(l;t)) + 

+ p(/(0;t) - /(i;0) + ?/ 2 (M) (22) 

From Eq.(22) we can now deduce the boundary condition for Eq.(21). Setting in Eq.(21) the 
correlation parameter A equal to 1 and comparing the terms in the rhs of Eq.(21) against 
the terms in the rhs of Eq.(12) we can infer that /(l; t) obeys 

l - (/(l;t) - /(0;t)) = P/(1;*) - ?/(2;«) (1 - /(!;*)) (23.o) 



Another pair of boundary conditions will be 

/(A;0) = C£P(X;t)p(X-\;t)) 

x 

and 



= P-, (23.6) 
t=o 



f(X^oo;t) - p_, (23.c) 

which mean that initially the lattice gas particles are uniformly distributed, with mean 
density p_, on the half-line X < 0, and that the density of the lattice gas at large separations 
from the phase boundary is equal to its unperturbed value. 

Equations (16) and (21) to (23) constitute a closed system of equations which allows a 
complete determination of X(t). Solution of these equations will be discussed in the next 
section. 
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IV. SOLUTION OF DYNAMICAL EQUATIONS IN THE CASE P+ = 0. 

We now turn to the continuous-space limit and rewrite our equations expanding /(A±l; t) 
into the Taylor series and retaining terms up to the second order in powers of the lattice 
spacing. We then find that f(X;t) obeys 

/(A;*) = \ ^/(A;t) - x(t) |^/(A;0, (24) 



while Eq.(23.a) transforms to 

= X(t)f(l;t), (25) 



kirx*™ 



A=l 

where, by virtue of Eq.(16), we have replaced the multiplier (p — q + qf(l; t)) by X(t). We 
note that Eqs.(24) and (25) hold for any relation between p and q (any orientation of the 
force F), but in the absence of the LDP the analysis of the case p > q does not make any 
sense. Consequently, in this section we will consider only the case when p < q. 

A. Expansion of the gas phase. 

Let us first consider the solution of Eqs.(24) and (25) supposing that X(t) > 0, (X(0) = 
0). Conditions when such a behavior takes place will be defined below. We notice that 
the structure of Eqs.(24) and (25) calls for the scaling solution in terms of variable uo = 
(A — 1)/X(t); < uj < oo. In terms of this variable Eqs.(24) and (25) can be rewritten as 

|1/h + (|a- W )(,-i)A /m = , (26) 

and 



>: 



= dx*(t)) f{u = 0), (27) 
„ at 



while Eqs.(23.b) and (23. c) collapse into a single equation 

f(oj = oo) = p_ (28) 
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Solution of Eqs.(26) to (28) can be readily obtained in an explicit form if we assume that 
dX 2 (t)/dt = A 2 , where A is a time-independent constant, < A < oo. Such an assumption 
actually makes sense if we recollect results of [f2^] and p0|-|32|, which demonstrated that in 



two extreme situations, i.e. when p/q = (totally directed walk of the PBP) and when 
p/q = 1 (no force exerted on the PBP), the PBP displacement shows the same generic 
behavior X{t) ~ \fi. Hence, one can expect that for arbitrary p/q, < p/q < 1, the PBP 
displacement should also grow in proportion to \ft. 

Time- independent A (p/q > 1). The general solution of Eq.(26) has the form 
/(w) = C x r dz exp(-^ (z 2 -2zj) + C 2 , (29) 
where C\ and C 2 are adjustable constants. Substitution of Eq.(29) into Eq.(27) gives 

d = A 2 C 2 , (30) 

while Eq.(28) yields the second relation 

f°° A 2 

d / dz exp(-— (z 2 - 2z)) + C 2 = p- (31) 
jo 2 

Consequently, we have for the density profile 

f{u) = ytj^a) {1 + A2 £ dz exp( ~T ^ 2 " 2 ^)>' ( 32 ) 

where the function I+(A) has been made explicit in Eq.(10). 

The density f(ui) in Eq.(32) is a function of A, which still remains undetermined. To 
define A we notice that X(t) ~ \ft behavior and Eq.(16) imply that 

f(u = 0) — >• , as t — > cxd, (33) 

and consequently, we have that in the limit t — >• oo the parameter A approaches a constant, 
time-independent value which obeys 

I + (A) = P-^-P^ (34) 
q-p 
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Equation (34) implicitly determines A as a function of p/q and p_. Numerical solution of 
this equation is presented in Fig. 3. 

Now, a simple analysis shows that Eq. (34) has a unique positive solution for any p and 
q which satisfy p > q(l — p_), (or, since p = 1 — q, such q which are less than 1/(2 — /?_)). 
When p/q — ► 1 — p_, the parameter A — ► as 

.4 „ ^ EZiilzi) , ,35) 

and is exactly equal to zero for p/q = 1 — p_ . It means that in the domain of parameters such 
that p > g(l — p_), the gas phase expands and the phase boundary moves as X(t) = Ayt, 
A > 0. Before we turn to the analysis of the behavior of the PBP in the regime p < q(l — p-), 
let us mention some other interesting aspects of Eqs.(34) and (32). 

Time-dependent A (p/q = 1). We note that Eq.(34) predicts that A diverges loga- 
rithmically when p/q — > 1 (Fig. 3). Nameley, 



A « J-2ln(l--), (36) 

which means apparently that when p = q the parameter A is some increasing function 
of time. This is, of course, consistent with the result of |E9l which states that the mean 



displacement of the PBP obeys X{t) ~ Jt ln(t) for the lattice gas with p_ = 1 and p = q = 
1/2. Let us now estimate, in terms of our approach, the behavior of A for arbitrary p_ and 
p = q = 1/2. For p = q our Eq.(16) reduces to 

X(t) = /(w = 0) (37) 

Next, supposing that X(t) still follows the law X(t) = A\fi, in which the prefactor A may 
be a slowly varying function of time, such that A/yi —*■ when t — >• oo, we find that the 
representation of /(A; t) in terms of a single scaled variable uj is still appropriate; weak time- 
dependence of parameter A actually results in the appearence of vanishing in time correction 
terms. We have then that the boundary condition in Eq.(25) reads 
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>: 



A 3 , 
« — ^=, when £ — > oo (38) 

a,=0 2 V* 



On the other hand, we can calculate the derivative of f{ui) directly, using Eq.(32). This 
gives 



>: 



P- A __, ,2 



.=0 V27T 



exp(-A 2 /2) (39) 



Comparing next the rhs of Eqs.(38) and (39), we infer that the parameter A obeys 



A 2 exp(A 2 /2) « J — , (40) 



which yields 



A « yin(^) (41) 



Equation (41) thus generalizes the result of [29] for arbitrary initial mean density p 



Wandering of the PBP in the critical case p/q = 1 — p _. Here we present some 
heuristic estimates of the time evolution of the second moment of the distribution P(X; t) 
in the case when the gas phase does not "wet" the region X > 0, i.e. when X(t) = 0. To do 
this, let us recall the Einstein relation between the diffusion coefficient D of a particle, which 
performs an unconstrained symmetric random walk in absence of external forces, and the 
mobility p of the same particle in the case when an external constant force is present. The 
Einstein relation states that p = (3D. Of course, it is not clear apriori whether the Einstein 
relation between the diffusion coefficient and the mobility should hold also for the tracer 
particle diffusing in a one-dimensional lattice gas; indeed, it may be invalidated because 



of the hard-core interactions. This question has been addressed for the first time in |33 
in which work several important advancements have been made. To illustrate some of the 
results obtained in |33|| , which are relevant to the model under study, let us first define the 



mobility of the tracer particle: 

p = lim t ^oop{t), (42) 



where //(£) denotes 



fi(t) = /m^o-^, (43) 



i.e. fi(t) is the ratio of the mean displacement X(t) of the tracer particle, diffusing in the 
presence of constant external force F, and F t; the ratio being taken in the limit when the 
external force tends to the critical value (zero) at which the mean displacement vanishes. 
Next, the diffusion coefficient of the tracer particle is defined by 



D = Um^ooD(t) = Umt->oc { X ^ glj) }, (44) 



where X^(F = 0,t) denotes the mean-square displacement of the tracer particle in the case 
when the external force is equal to its critical value (zero). Now, for the tracer diffusion in 



a one- dimensional hard-core lattice gas one has that X^(F = 0, t) obeys Eq.(l), while X(t) 
is determined by Eqs.(3) and (4). 

One readily notices now that the Einstein relation holds trivially for infinitely large 



systems, since here both p, and D are equal to zero [33|. A more striking result obtained in 
PB| concerned the case when the one- dimensional lattice is a closed ring of length L. It was 
shown that here both p and D are finite, both vanish with the length of the ring as 1/L 
and obey the Einstein relation fi(L) = (3D(L) exactly! Next, |3l| and subsequently, []32 



focused on the non-stationary behavior in infinite systems and showed that the Einstein 
relation holds in an even more general sense: namely, the time-dependent mobility //(£) and 
the diffusivity D(t) obey 

M*) = PD(t), (45) 

at times t sufficiently large, such that the asymptotical regimes described by Eqs.(l) and 
(3) are established. 

Now, in the situation under study we have non-zero critical force (A = when p = 
q(l — p-) or, in other words, when F = F c = — fi~Hn{\ — p_)). We thus define the time- 
dependent mobility p(t) as 
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X(t) 
fi(t) = lim F ^ Fc ——j—, (46) 

which yields, by virtue of Eq.(35), the following result 

"W = " ^ {I (47) 

Assuming next that the generalized Einstein relation in Eq.(45) holds in the situation under 
study, we find that in the critical case p — q(l — p_) the mean-square displacement of the 
PBP obeys: 

I- p. 



X*(F = F e ,t) = -\~i (4f 

P- V 7T 

which is surprisingly similar to the classic result in Eq.(l). 

In Fig. 4 we compare our analytical prediction in Eq.(48) against the results of Monte 
Carlo simulations, performed at three different values of the gas phase densities p_. It 
shows that our Eq.(48) is in a good agreement with the numerical results. This means that 
the Einstein relation in Eq.(45) holds even in such a "pathological" situation, in which the 
critical value of the external force is not equal to zero and the particle density is different 
from both sides of the test particle. 

Density profiles. Let us now analyse the form of the density profiles as seen from the 
PBP. In Figs. 5 and 6 we plot the result in Eq.(32) versus the scaled variable u for different 
initial mean densities and different values of the ratio p/q. 

In Fig. 5 we depict f(oo) for fixed p/q = 0.9, which corresponds to fixed "boundary 
tension" force F, and different initial mean densities p_. In the range of used parameters, 
all corresponding values of A are of the same order (A rj 1) and the density curves look 
quite similar; starting from the same value at uj — 0, f(cu = 0) = 1 —p/q = 0.1, they 
quite rapidly, within a few units of uj, approach their unperturbed initial values. We note, 
however, that on the X-scale it does not mean that the density past the rightmost particle 
rapidly reaches the unperturbed value p_. Instead, at sufficiently large times the density 
stays almost constant and equal to 1 —p/q within a macroscopically large region ~ X(t) . 
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Now, in Fig.6 we plot f(u) versus uo in the opposite case when p_ is fixed and the 
"boundary tension" force is varied. Here the density profiles display rather strong depen- 
dence on the parameter A. When A is small f(u) shows almost linear dependence on uo 
(curves (3) and (4)). The reason for such a behavior is that here the phase boundary moves 
essentially slower (A < 1), compared to the typical displacements of the gas particles, which 
then have sufficient time to equilibrate the density profile past the PBP. In the opposite case 
of relatively large values of A, (A > 1), such an equilibration does not take place and the 
dependence of f(oo) on uo is progressively more pronounced the larger A is. 

It may also be worth-while to discuss the shapes of the density profiles in terms of the 
variables A and t. First, from Eq.(32) we have that in the limit of small oo, i.e. A <C X(t), 
the density obeys 

f(X;t) «(l-£)[l+ A (A ~ 1} + ...], (49) 
q V* 

which means that past the PBP the density is almost constant in the region whose size 
grows in proportion to X(t). Next, within the opposite limit, i.e. at distances A which 
exceed considerably X(t), we obtain from Eq.(32) the following result 

f(X,t) « p _ - (l-£) ^ exp(-^) + ... (50) 

Equation (50) shows that at large separations from the phase boundary the density ap- 
proaches the unperturbed value p_ exponentially fast. The approach is from below and is 
weakly (only through the prefactors) dependent on the parameters p/q and A. 

Mass of particles and mean density. We close this subsection with a brief analysis 
of the time evolution of the integral characteristic of the propagating gas phase; namely, of 
the "mass" M(t) and the mean density p mean = M(t)/X(t) of lattice gas particles at sites 
X > at time t. 

The parameter M(t), which measures the amount of the gas-phase particles which 
emerged up to time t in the previously empty half-line X > 0, is formally defined as 
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rXit) 

M(t) = / dX p(X;t) (51) 

JO 



Changing the variable of integration, we find that M(t) can be rewritten as 

M(t) = / d\f(X;t) = 

Jo 



where M is given by 



X(t) f dw fho) = Mt 1/2 , (52) 
Jo 



M = A(l-Z) exp(A 2 /2) (53) 
Q 



Figure 7 displays the plot of the prefactor M versus p/q and shows that M is a mono- 
tonically increasing function of p/q. In contrast to the parameter A, M remains finite for 
p = q, which means that bulk contribution to the "mass", as it could be expected intuitively, 
comes from the lattice gas particles, whose motion is constrained by hard-core exclusions 
and whose mean displacement grows only as \/t, without an additional logarithmic factor 
which is specific only to the PBP. 

Finally, we depict in Fig.7 the mean density on the interval X e [0, X(t)] , defined as 

Pmean = ^ = (1 " l ) CX P (A 2 /2) (54) 

X[t) q 

Figure 7 shows that despite of the exponential factor exp(A 2 /2) the mean density p mean 
rapidly decreases with an increase of p/q and is a slowly increasing function of p_. 

B. Compression of the gas phase. 

Let us next address the question of the PBP dynamics in the case p < q(l — p_), when 
X{t) is expected to be less than zero and thus the gas phase to be effectively compressed 
by the "boundary tension" force exerted on the PBP. Recollecting the results of [|3(J-|32| 
we suppose that here X(t) obeys X(t) = —By/i, B > 0, and define the scaled variable as 
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uo = (A — l)/By/i, where uo is positive definite < uo < oo. In terms of this variable Eqs.(24) 
takes the form 

d 2 ._ - , ^ d 



f(u) + B 2 (uo+l) — f(uo) = 0, (55) 



duo 2 duo' 

while the boundary condition in Eq.(25) reads 



>■ 



B 2 f(uo = 0) (56) 



ui=Q 

Again, the boundary and initial conditions in Eqs.(23.b) and (23. c) collapse into a single 
Eq.(28). 

The general solution of Eq.(55) can be written down as 

/M = C, r dz exp(-^(z 2 + 2z)) +C 2 , (57) 
jo 2 

where C\ and C 2 are to be chosen in such a way that Eqs.(28) and (56) are satisfied. Inserting 
Eq.(55) into Eqs.(28) and (56) we then obtain 

d = - B 2 C 2 , (58) 

and 

/■oo R 2 



C 2 = p_{l - B 2 / dz exp(- — (z 2 + 2z))}- 1 (59) 

Consequently, the density profile past the PBP can be expressed in terms of B and uo as 

f{uo) = p-{\- B 2 r dz exp(-^V + 2z))}/ 

Jo 2 

/•oo R2 

/{l - fi 2 / dz exp{- — {z 2 + 2z))} (60) 
jo 2 

Next, Eq.(16) implies that also in this case f(uo = 0) — ► (q — p)/q as t — > oo, which yields 
eventually the following closed-form equation for the parameter 5: 



-B 2 / 9 u QP- 

o 2 g — p 



I - B 2 / dz exp(- — (z 2 + 2z)) = -^=- (61) 
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Equation (61) can be put into a more compact form if we express the integral over dz in 
terms of the probability integral. We then obtain 

J_(B) = ^-P-^-P , (62) 
q-p 

where I -(B) is defined in Eq.(lO). In Fig.8 we present the numerical solution of Eq.(62), 
plotting the prefactor B as a function of the ratio p/q at different values of the density p_. 

Equation (62) resembles the form of Eq.(34), which determines the parameter A, but 
differs from it in two aspects; first, the rhs of Eq.(62) is exactly the rhs of Eq.(34) but taken 
with the opposite sign, which insures that B is positive for p < q(l —/?-), and second, the 
sign before the probability integral in brackets is opposite to that in Eq.(34). The latter 
circumstance is responsible for the fact that B tends to the limiting value Bu m when q — > 1 
(p — > 0). When q — > p/(l — pJ) the parameter B tends to zero exactly in the same way as 
the parameter A in Eq.(35) taken with the opposite sign, which means that the prefactor 
in X(t) does not have a discontinuity at the "critical" point p/q = 1 — p_ both for its value 
and for its slope. 

Now, in the limit q = 1 (p = 0), Eq.(62) reduces to 



2 B Um exp(^ m /2) [1 - $(Aim/V2)] = l-p_, (63) 



which was obtained previously in [p0l , |3Tfl (see Eq.(2) in the present paper). Within the limit 
p_ —> 1 Eq.(63) yields 

B lim « y^(l-p-), (64) 

which shows that Bi im , as it could be expected intuitively, tends to zero when the density 
tends to 1. 

When the gas is very dilute, i.e. p_ <C 1, we may expect that Bi im is large. Expanding 
the probability integral as 

®(B lim /V2) « 1 - y ^ Bf im exp(-B? im /2) + 

24 



-B H l exp(-5L/2) - .- , (65) 

7T 

we find, upon substitution of Eq.(65) into the Eq.(63), the following result 

Bum w -— , (66) 

i.e. Bii m diverges when p_ — > in proportion to the inverse of the square-root of the particle 
mean density. In Fig. 9 we present the numerical solution of Eq.(63) together with the results 
of Monte Carlo simulations. Obviously, the agreement is very good. 

Finally, in Fig. 10 we combine the results of the subsections A and B and plot both analyt- 
ical and Monte Carlo results obtained for the dependence of the prefactor a(F) = X(t)/y/i 
on the ratio p/q and the density p_. Again, we find very good agreement between our ana- 
lytical predictions and numerical results, which support the validity of the approximations 
involved in our analysis. 

Density profiles. Consider now the density profiles as seen from the PBP in the 
compression regime. In Fig. 11 we plot f(u) versus u> for different values of p/q at fixed p_. 

Figure 11 shows that similarly to the behavior in the expansion regime, the density 
profiles are quite sensitive to the value of the parameter B. When B is smaller than unity, 
f(ui) shows almost linear dependence on u, while in the case when B > 1 this dependence 
is non-linear and f(u) rapidly drops from f(u = 0) = 1 — p/q to p_. 

We finally present explicit results for f(X;t). In the limit of small A, such that A <C 
X(t)/B 2 , we have 

/(A;*) »(1-5)[1- ^^i> +...], (67) 

q yt 

which shows that the density is almost constant, (being only slightly less than (1 —p/q)), 
in the spatial region whose size is of the order of the PBP mean displacement. 
For large A we obtain from Eq.(60) 

n R\/2t 
f(\;t) » p_ + (1 - q -^— exp(-\ 2 /2t) + ... , (68) 

q A 
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i.e. similarly to the behavior in the expansion regime, the density approaches the un- 
perturbed value p_ exponentially fast and with a rate which is weakly (only through the 
pre-exponential factor) dependent on the parameter B and the ratio p/q. 

V. DYNAMICAL EQUATIONS IN THE PRESENCE OF THE LOW-DENSITY 

PHASE. 

Let us now consider the time evolution of the local density p(X; t) and of the probability 
distribution P(X; t) in the general case when the LDP is present and p_ > p + > 0. 

One readily notices that also in this case Eq.(12) and, consequently, Eq.(13), describe 
the time evolution of the realization-averaged occupation variable p(X; t) for all X excluding 
the sites X = X(t) ± 1. Dynamical equation describing evolution of p(X; t) at X = X(t) — 1 
will be, however, somewhat modified as compared to Eq.(14). We have here 

p(X(t)-l;t) = l -{p{X{t)-2 ] t) - p (X(t)-l;t)) + 
+ q(l - p(X(t)-l;t))p(X(t)-2;t)- 

- p (1 - p(X(t) + 1; t)) p(X(t) - 1; t), (69) 

in which we account that the hops of the PBP in positive direction can be constrained by 
the LDP particles by introducing a factor (1 — p(X(t) + 1; £)). In a similar fashion, we find 
that at the site X — X{t) + 1 the local particle density obeys 

p(X(t) + l;t) = l -{p{X{t) + 2 ] t) - p(X(t) + l;t))- 
-q(l - p(X(t)-l;t))p(X(t) + l;t) + 

+ p (1 - p(X(t) + 1; £)) p(X(t) + 2; t) (70) 

Next, for the time evolution of the distribution function P(X; t) we obtain the following 
equation 

26 



P(X;t) = -P(X;t)\p(l - p(X + l;t)) + q (1 - p(X - l;t))} + 

+ (1 - p(X; t)) [p P(X - 1; t) + q P(X + 1; £)], (71) 

which differs from the corresponding equation of the previous sections, Eq.(15), by the 
factors (1 — p(X + 1;£)) and (1 — p(X;t)) in the first and third terms respectively; these 
factors account, in a mean- field-type fashion, for the fact that hops of the PBP in the positive 
direction can take place only if the corresponding lattice sites are free of the LDP particles 
at this moment of time. 

Further on, multiplying both sides of Eq.(71) by X and summing over all lattice sites we 
have that the mean displacement of the PBP obeys: 

X (t) = p - q - pf(\ = -l;t) + qf(\ = l;t), (72) 

which thus generalizes Eq.(16) for the case of non-zero density of the LDP; the factor /(— 1; t) 
on the right-hand-side of Eq.(72) accounts for the hindering effects of the LDP particles on 
the PBP dynamics. 

Consider now the time evolution of the correlation function /(A;t), defined in Eq.(17). 
By virtue of Eqs.(18), (13) and (71) we find that the evolution of this property is guided by: 

/(A;i) = i[/(A + l;t) + f(X-l;t) - 2 f(X;t)] - 
-f(X;t)[l - pf(-l;t) - qf(l;t)} + 

+ p /(A - 1; t) [1 - /(-l; £)] + q /(A + 1; t) [1 - /(l; £)], (73) 

which holds for all A excluding A = ±1. In the limit p_ — ► 0, i.e when /(— l;i) — ■* 0, this 
equation reduces to Eq.(21). In the continuous-space limit Eq.(73) attains the form 

/(A;*) = l^ 11 - ip-Q-pf(-l',t)+ q ffrt)]°i^ (74) 

which is exactly Eq.(24) with X(t) defined by Eq.(72). 
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Further on, we find that the correlation function /(A; t) at the left-hand adjacent to the 
PBP site (for A = 1) obeys 



/(!;£) = - (/(2;t) - /(!;«)) + p/(0;t) (1 - /(-!;*)) - 



g/(l;t)(l - /(!;*)) - 2p/(l;i)(l - /(-!;*)) + 



+ 2g/(2;f)(l - /(l;i)) (75) 
Comparing now Eq.(75) with Eq.(73) we have the following condition on f(X;t) at A = 1: 

l - (f(0;t) - /(1;0) = qf(2;t) (1 - f(l;t)) - p f(l;t) (1 - /(-l;t)) (76) 
Next, from Eqs.(18), (70) and (71) we can derive 



f(-l;t) = -(/(-2;t) - /(-!;*)) 



-p/(-l;i)(l - /(-!;*)) - 2g/(-l;i)(l - /(!;*)) + 



+ 2pf(-2;t)(l ~ /(-!;*)) + 9/(0;*)(l " /(I;*)), (77) 



which allows us to deduce the boundary condition on /(A; t) at the point A = — 1: 



-(/(0;t) - /(-!;«)) = -g/(-l;*)(l " /(!;«)) + 



+ p/(-2;t)(l - /(-!;*)) (78) 



In the continuous-space limit Eqs.(76) and (78) reduce to 

1 df(X;t) 



2 d\ 



X(t) f(±l;t), (79) 



A=±l 



which represent two boundary conditions for the continuous-space Eq.(74). Equations (74) 
and (79), with the initial conditions 
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/(A;*)| t=0 = P+forA<0, (8O.0) 

/(A;t)lt=o = P-forA>0, (80.6) 
and the boundary conditions 

/(^)Ia^oo = P-, (81-a) 

/(A;*)Ia— 00 = P+. (81-6) 

constitute a closed system of equations which allows to compute X(t) and the density profiles 
for arbitrary relation between p and q, as well as for arbitrary p + and p_. 

VI. SOLUTION OF DYNAMICAL EQUATIONS IN THE GENERAL CASE 

P- > P+ > 0. 

In this section we will derive explicit results for the dynamics of the mean displacement 
of the PBP and also for the density distribution around it. As it was done in the previous 
sections, we will discuss separately the behavior in the case when the HDP expands, com- 
pressing the LDP, and when, on the contrary, the LDP and the external force F compress 
the HDP. 

A. Expansion of the high-density phase. 

We again set X(t) = Ay/i and suppose first that A > 0. Conditions at which such a 
behavior takes place will be specified below. For A > 1 (past the PBP) we then have 

A 2 

f(u) = 1 + ^ {1 + A 2 £dz exp(- T (z 2 -2z)}, (82) 

where u = (X—l)/A^/i and I+(A) is defined in Eq.(10). In front of the PBP, i.e. for A < — 1, 
the scaled density profile is given by 
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f{6) = 1 _ P I + {A) (I - # [dz exp(-^(z 2 + 2z)}, (83) 

in which we have denoted 6 = — (A + l)/A\/i and I -(A) is made explicit in Eq.(lO). The 
density distributions f(u) and f(6) for different values of the parameters A, p± and p (q) 
are depicted in Figs. 5, 6 and 11 respectively. 

Equations (82) and (83) contain the parameter A, which has not yet been specified. To 
determine A we take advantage of Eq.(72) which yields the following condition on the local 
densities at the sites adjacent to the PBP position: 

9(1 - /(A = l;t)) = p(l - /(A = -l;t)) (84) 

Upon substitution of Eqs.(82) and (83) into the latter equation we find that A (in case when 
A > 0) obeys the following transcendental equation: 

«P- PP+ - a - v (85) 

1 + I + (A) l-I-(A) ~ q V " [Hb) 

which generalizes our Eq.(34) and also the result of |H] (Eq.(4) of the present paper) for 
the case when the particle densities from the left and from the right of the PBP are different 
and the density of the LDP is not zero. One directly verifies that Eq.(85) reduces to Eq.(34) 
when we set p + = 0, while setting p + = p_ we recover Eq.(4). 

Let us now find the conditions under which the parameter A is positive, i.e. the HDP 
expands. To do this, we simply notice that when A = both I+(A) and I -(A) are equal to 
zero, which means that the "critical" relation between p, q and p± is: 

q (1 - p_) = V (1 - P+) (86) 

Equation (86) implies that A vanishes, (i.e. the LDP and the HDP are in equilibrium with 
each other), when the probability of the PBP to go towards the HDP times the density of 
vacancies in this phase is exactly equal to the probability of going towards the LDP times 
the density of vacancies in this phase. When p(l — p + ) > q(l — p_) the HDP expands. 



We note that Eq.(86) was previously obtained in |33[ from the analysis of the stationary 
states in a one-dimensional lattice gas placed in a finite box of length L. By explicit calcula- 
tion of the distribution function of the PBP position in the general case when the numbers 
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of the lattice gas particles from the right and from the left of the PBP are not equal, it was 



found |33| that in the limit L — > oo the PBP is localized at point X , which divides the 
system in proportion given by Eq.(86). 

Equation (86) can also be rewritten using the definition of the external force F. Upon 
some algebra, we find then that the critical force at which both phases are in equilibrium 
with each other is given by 

F c = /T 1 ln{\^±) (87) 
1-p- 

Now, let us discuss the behavior of the parameter A in the limit when A is small or large, 
and calculate the diffusivity of the PBP in the critical case A = 0. In the limit of small A, 
i.e. when p, q and p± are close to their "critical" values determined by Eqs.(86) and (87), 



both I±(A) ~ Jtt/2A. Substituting these expressions into Eq.(85) we find 

V vr qp-+ pp+ 

which is valid when A<1. Eq.(88) allows for the computation of the PBP mobility, which 
we determine following the arguments presented in Section IV as 

X(t) 



p(t) = lim 



F^F C 



(F - F c ) t 



= r 1 ' 2 Um F ^ Fc {F _ Fc) (89) 
Substituting Eq.(88) into Eq.(89) and taking the limit F — > F c , we find 

(p- +P+- 2p-P+) V ?rt 
which yields, by virtue of Eq.(45), the result presented in Eq.(7). Eq.(7) generalizes the 
classical result in Eq.(l) for the situation in which the mean particles densities for both 
sides of the tracer particle are different from each other. One can directly verify that Eq.(7) 
reduces to Eq.(l) when p_ = p + , while setting p + = we recover our previous result 
in Eq.(48). In Fig. 12 we compare our analytical prediction in Eq.(90) against the results 
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of Monte Carlo simulations, which shows that an approximate approach developed here 
represents a fair description of the PBP dynamics. 

Next, in Section IV we have demonstrated that the prefactor A diverges when p — > q. 
Consequently, we can expect that even in the presence of the LDP the prefactor A can attain 
large values when p + <C 1 and p/q — ► 1. Setting in Eq.(85) q = p and using the expansion 
in Eq.(65) we find from Eq.(85) that A is defined in the limit p + — ► by 



A » y/2ln(p-/p + ), (91) 



i.e. A grows as a square-root of the logarithm of p + when p + — ► 0. 

Finally, we estimate the behavior of the ratio <5 of the particle densities immediately past 
and in front of the PBP. At zero moment of time this ratio is evidently S = Sq = p-/p + . 
Our results in Eqs.(82) and (83) suggest that after some transient period of time the density 
profiles past and in front of the PBP attain stationary forms with respect to the variable 
uj = (A — 1)/ X(t). Consequently, we have that, as the time evolves, the ratio of the particle 
densities immediately past and in front of the PBP tends to a constant value 

1-I-(A) 

6 = 6 ° TTJ^A) (92) 

Eq.(92) holds for arbitrary values of A. In the asymptotic limits when A is small or large, 
we find from Eq.(92) the following explicit asymptotic forms for 5: 

5 « 5 (1 - (- - 1) A 2 + ... ), when A < 1, (93) 



and 



2 



^ exp(A 2 /2) 

V2TA 3 V ; 



Therefore, the parameter 8, as it could be expected intuitively, is always less than Sq. Com- 
plete dependence of S on the parameter A is presented in Fig. 14. 
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B. Compression of the high-density phase. 

Consider now the behavior in the regime when the LDP and the applied force compress 
the HDP. Setting X(t) = —B\/i, where B is supposed to be a positive constant, we find 
from Eqs.(74) and (79) that the particle density for A > (i.e. at sites X < X(t)) obeys 
Eq.(60), in which the variable to is defined as 10 = (A — l)/By/i and the parameter A is 
replaced by B. From the other side of the PBP, i.e for A < 0, we have 

/(*) = i + P j +{B) (I + B 2 J°dz exp(-|V-2^))}, (95) 

where the scaled variable 6 = — (A + l)/By/i. The function f{6) is depicted in Figs. 5 and 
6. Substituting Eqs.(60) and (95) into Eq.(84) we arrive at the following transcendental 
equation for the parameter B: 

qP ~ PP+ -Q-P 06) 



1 - I_(B) 1 + I + (B) 

Eq.(96) thus generalizes the result in Eq.(62) for the case of the non-zero particles density 
in the LDP. 

Now, noticing that Eqs.(96) and (85) can be cast into one another by the substitution 
±I±(A) — ► ^fI T (B), we can construct a general equation for the parameter a(F) in Eq.(3). 
This equation is presented in Eq.(6) and holds for arbitrary relation between p± and p/q, 
describing hence both the expansion and the compression regimes. In Fig. 13 we present the 
numerical solution of Eq.(6), plotting a(F) as a function of the ratio p/q for different values 
of p_ and p + . 

Finally, we analyze the behavior of the parameter S, which is defined as the ratio of the 
particle density immediately past the PBP and the particle density immediately in front of 
the PBP. Using Eqs.(60) and (95) we find 

! + /+(£) 

6 = 6o T^rm (97) 

Numerical plot of S(B) is presented in Fig. 14. 
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Asymptotic behavior of the parameter 5 in the limits when B is small or large readily 
follows from our Eqs.(93) and (94). Here we have 

5 « 5 (1 + (| - 1) B 2 + ... ), when S < 1, (98) 

and 

5 ^ 5 V^ B 3 exp(B 2 /2), when B > 1, (99) 

which means that in the compression regime the parameter <5 is always greater than Sq. 

VII. MONTE CARLO SIMULATIONS. 

In order to check our analytical predictions, derived in terms of a mean-field approxima- 
tion, we have performed Monte Carlo simulations of the process defined in the begining of 
Section II. The simulation algorithm was defined as follows: 

We constructed first a one-dimensional regular lattice of unit spacing and length 2L + 1, 
sites of which were labelled by integers of the interval [—L,L]. In all simulations we took 
L = 10 3 . At the zero moment of the MC time the particles were placed randomly on the 
lattice with the prescribed mean densities and the constraint that two particles can never 
simultaneously occupy the same site. To do this, we have called, for each lattice site from 
the interval [-L + 1, —1] independently, a random number from the interval [0, 1]. In case 
when the random number produced by the generator was less that p_ a particle was created 
on this lattice site. In case when the random number was greater than p_ the site was left 
empty. The same routine was performed for the sites with positive numbers [1, L — 1]; here 
a particle was created at the corresponding site in case when the random number was less 
than p + and the site was left empty if the random number was greater than p + . The phase 
boundary particle was placed at the origin. Additionally, we have prescribed that the sites 
X = ±L are occupied by particles. The particles at theses sites X = ±L are made immobile, 
blocking the lattice from both sides and preventing other particle to leave the system. 
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The subsequent particle dynamics employed in our simulations follows the definitions of 
Section 2 closely. We call for a random integer number X from the interval [— L + 1, L — 1]. 
Here three different events may take place: 

(i) If the site X is vacant, a new site is considered. 

(ii) If the site X is occupied by a particle, we first increase the MC time by unity and 
then let the particle choose, at random, a potential jump direction. This is done again by 
calling a random number from the interval [0, 1]. If the random number is less than 0.5, 
the particle attempts to jump to the site X — 1; otherwise, it attempts to jump to the site 
X + 1. The jump is fulfilled if at this moment of the MC time the adjacent site in the chosen 
direction is vacant (not occupied by any other particle or the PBP). Otherwise, the particle 
remains at X . 

(iii) If the site X appears to be occupied by the PBP, we increase the MC time by unity 
and consider a random number from the interval [0, 1]; in case when this number is less than 
the prescribed value q, the PBP attempts to jump to the site X — 1. Otherwise, it attempts 
to jump to the site X + 1. The jump is fulfilled if at this moment of MC time the adjacent 
site in the chosen direction is vacant. Otherwise, the PBP remains at X. 

In simulations we have followed the time evolution of several different properties: the 
PBP displacement, squared displacement (under the critical conditions) and the occupations 
of the sites X = ±L =F 100. Time behavior of these properties was plotted versus the 
"physical time" t, which is the time needed for each particle to move once, on average, or 
in other words, t = MCtime/number of particles. We have observed that for all values of 
the parameters p± and q, used in our simulations, the stationary regime in which the ratio 
X r {t)/\/t approaches a constant value is established for displacements of order of 200 lattice 
units. To get the spatially resolved behavior, in computation of the PBP displacement 
each realization of the process was interrupted at the moment when the absolute value 
of the PBP displacement reaches the value of 500 lattice. For calculation of the mean 
displacement we used typically 10 2 realizations for each set of parameters p± and q. Results 
of these simulations are presented in Figs. 9 and 10. Further on, computing the squared 
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displacement we interrupted each realization of the process at the moment when the span 
of the PBP trajectory is equal to 10 2 . Mean-square displacement was obtained by averaging 
over 2 x 10 3 realizations. Results for mean-square displacement of the PBP are presented 
in Figs. 4 and 12. In all cases, we have obtained remarkably good agreement between our 
analytical predictions and simulation results. Finally, the measurements of the occupations 
of the sites X = ±L =F 100 were performed in order to be sure that the perturbances created 
by the PBP do not spread during the simulation time through the whole system and do 
not lead to artificial behaviors associated with the finite-size effects. We have observed 
that actually the mean densities of these sites don't vary with time and are equal to the 
unperturbed values p±. 

VIII. CONCLUSIONS. 

To conclude, we have examined in terms of a mean-field- type approach the dynamics 
of the phase boundary propagation in a one-dimensional hard-core lattice gas which was 
initially put into a non-equlibrium, "shock"-like configuration and then allowed to evolve 
in time by particles attempting to hop to neighboring unoccupied sites. The "shock" con- 
figuration means that particle mean densities from the left and from the right of the origin 
are different. All particles of the lattice gas, except the particle separating the high- and 
the low-density phases, have symmetric hopping probabilities, while the phase boundary 
particle is subject to a constant force F and has asymmetric hopping probabilities. We have 
shown that the mean displacement of the PBP follows the generic law X{t) = a(F)\/i, in 
which the parameter a(F) can be both positive and negative, depending on the relation 
between the magnitude of the force and the initial mean densities. This prefactor is deter- 
mined implicitly, in a form of the transcendental Eq.(6) for arbitrary magnitude of the force 
and arbitrary relation between the particle densities in the high- and low-density phases. In 
several asymptotic limits we find explicit formulae for the prefactor. Further on, we have 
shown that when F is equal to the critical value F c , Eq.(87), the parameter a(F) is exactly 
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equal to zero. In this case the high- and the low-density phases are in equilibrium with each 



other. We have found that here the mean-square displacement X%(t) of the PBP follows 



X^(t) ~ r )\fi, i.e. shows a sub-diffusive behavior. The form of the prefactor 7 is determined 
explicitly, Eq.(7). Our analytical findings are in a very good agreement with the results of 
numerical simulations. 

The authors wish to thank J.L.Lebowitz and R.Kotecky for helpful and encouraging dis- 
cussions. Financial support from the FNRS and the COST Project D5/0003/95 is gratefully 
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Fig. 1. Initial "shock" configuration of the lattice gas particles: p_ is the mean density of 
particles at the half-line — oo < X < and p+, (p+ < p_), is the mean density of particles 
at the half-line < X < oo. All gas particles, excluding the "phase boundary" particle, have 
equal (= 1/2) probabilities for jumps to the right and to the left. The PBP has asymmetric 
hopping probabilities. 




Fig. 2. Initial "shock" configuration of the lattice gas particles. All particles are placed at 
random positions and at a fixed mean density p— from the left to the origin of the lattice, 
i.e. at sites X < 0. All sites X > are vacant at t = 0. 




Fig. 3. Numerical solution of Eq.(34). The parameter A is plotted versus the ratio p/q for 
different values of the mean density p_ . The curves from the left to the right correspond to 
p- = 0.9,0.8,0.7, ... ,0.1. 
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Fig. 4. Mean-square displacement of the PBP in the critical case p/q = 1 — p- . The solid 
lines depict our analytical prediction in Eq.(48). The noisy lines give the corresponding 
results of Monte Carlo simulations. The curves from the top to the bottom correspond to 
P- = 0.3, 0.5 and 0.7 respectively. 




Fig. 5. Expansion of the gas phase. Plot of /(<*>), Eq.(32), versus ui — (A — 1)/X(t) at a 
fixed ratio p/q and variable p_ . The corresponding values of the parameter A for the curves 
from the top (p_ = 0.9) to the bottom (p_ = 0.5) are A = 1.37, 1.3, 1.25 and 1. 
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Fig. 6. Expansion of the gas phase. Plot of f(ui), Eq.(32), versus ui at fixed density p_ = 0.9 
and variable ratio p/q. The corresponding values of the parameter A are 0.09 (4), 0.31 (3), 
0.58 (2) and 1 (1). 




Fig. 7. Expansion of the gas phase. Plots of the parameter M and the mean density p m ean 
in the half-line X > versus pjq for different initial densities p_ . The curves (a) to (e) show 
the value of the parameter M for p_ = 0.9,0.8,0.7,0.6 and 0.5 respectively. The curves (1) 
to (5) display the corresponding mean densities p m ean- 




Fig. 8. Compression of the gas phase. Numerical solution of Eq.(62). The parameter B is 
plotted versus the ratio pjq for different values of the mean density p_ . The curves from top 
to bottom correspond to p_ = 0.1,0.2,0.3, ... ,0.9. 




Fig. .. Plot of the parameter B n m versus the mean particle density p- . The solid line 
shows the numerical solution of Eq.(63). The diamonds represent the associated Monte 
Carlo simulation results. 




Fig. 1.. Theoretical and experimental results for the prefactor in the dependence X[i) — 
a(F)^/t as the function of the ratio p/q. The curves from top to bottom correspond to 
different values of the density p—\ the upper curve gives a(F) for p— — 0.6, the lower curve 
corresponds to p- = 0.4 and the curve in the middle - to p- = 0.5. Symbols denote the 
results of Monte Carlo simulations. 
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Fig. 11. Compression of the gas phase. Plot of /(<*>), Eq.(60), versus the scaled variable w at 
fixed density p_ = 0.2 and variable ratio p/q. The corresponding values of the parameter B 
are 1.52 (1), 1.21 (2), 0.83 (3) and 0.34 (4). 




Fig. 12. Mean-square displacement of the PBP in the critical case p/q = (1 — /a)/(l — P+), 
Eq.(86). The solid lines show our analytical prediction from Eq.(7) and the noisy lines give 
the results of Monte Carlo simulations. The curves from top to bottom correspond to the 
following values of the parameters (p_ ,p+): The first two curves are the analytical result in 
Eq.(l) and numerical data for the symmetric cases (0.4,0.4) and (0.5,0.5). The lower curves 
correspond to (0.6,0.5), (0.9, 0.4), (0.7, 0.5) and (0.8,0.2) respectively. 
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Fig. 13. Numerical solution of Eq.(6). Prefactor a(F), Eq.(3), is plotted versus the ratio 
p/q for different initial mean densities from the left and from the right of the origin. The 
curves from (1) to (4) correspond respectively to (p- = 0.9, p+ - 0.3), (0.8,0.2), (0.7,0.3) 
and (0.5,0.5). 




Fig. 14. Plot of the parameter 5, Eqs.(92) and (97), versus the parameters A and B. The 
curve (1) describes the behavior of 5(B) in the compression regime. The curve (2) defines 
5(A) in the expansion regime. 



